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Abstract. We use nonlinear kinetic theory to study the remnant dynamics and the particle acceleration as well as the properties 
of the nonthermal emission from the supernova remnant SN 1006. The known range of astronomical parameters is examined to 
determine whether it encompasses the existing synchrotron emission data. Given the present-day spatial extent and expansion 
rate of the object, it is shown that the hadronic gamma-ray flux is very sensitive to the ambient gas density Nu and that the 
existing H.E.S.S. upper limit requires A'h < 0.1 cm"'. The strength of the amplified magnetic field downstream of the shock is 
about 150 yuG. 
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1. Introduction 

I Despite the widespread belief that the cosmic rays (CRs) below 
' the so-called knee in the energy spectrum at a few 10'^ eV are 
_ accelerated in the shell-type supernova remnants (SNRs) of our 
■ Galaxy, this proposition still has only limited evidence. This is 
' not due to an inadequate theoretical understanding. In fact, es- 
sentially all physical processes involved have either been clar- 
ified or, Uke the injection rates of ions and of electrons into 
the diffusive shock acceleration process, can be inferred from 
the observed radio and X -ray synchrotron spectra (for a recent 
overview, see IVoIkl2004l) . The Umitation is rather the scarcity 
of SNR detections in TeV y-rays as an unequivocal signature 
of the existence of very high energy particles in these sources. 

Of this SN 1006 (G327.6-H14.6) is the most promi- 
nent case. Observations of nonthermal X-rays feovam a et alJ 
Il99 5) suggest that at least CR electrons are accelerated in 
SN 1006 up to energies of about 100 TeV. Su bsequent y-ray 
observations with the CANGAROO telescopes JTanimori et al.l 
[ 1998, 2001) strengthe ned this conclusion. It was also shown 
Iserezhko et alJl2002l) that nonlinear kinetic theory of CR ac- 
celeration is consistent with all observational data for a value 
A^H = 0.3 cm"^ of the ambient interstellar medium (ISM) den- 
sity from the range 0.05 < A^h ^ 0.3 cm"^ existing in the 
literature. However, SN 1006 could not be detected by the 
H.E.S.S. experiment as a TeV source in a total of 18. 2h (in 
2003) and 6.3h (in 2004) livetime of ON source observations 
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with SN 1006 in the field of view jAharonian et alJ2005l) . The 
H.E.S.S. upper limit is roughly one order of magnitude lower 
than the published CANGAROO flux. 

In this paper we demonstrate that the H.E.S.S. upper 
limit does not invalidate the theoretical picture jBerezhko et alJ 
I2OO2) on which the previous calculation of the y-ray emission 
spectrum has been based. As we shall show, it is rather the 
value of the external astronomical parameter A^h that strongly 
influences the hadronic y-ray flux. The Inverse Compton (IC) 
y-ray flux, on the other hand, depends only weakly on A^h- We 
use this theory to describe all the relevant properties of SN 1006 
and re-examine the most relevant set of physical and astronom- 
ical parameters, mainly A^h, the total explosion energy and, to a 
minor extent, the interior magnetic field strength as well as the 
distance to the object, given the H.E.S.S. upper limit. 

2. SNR parameters 

Since SN 1006 is a type la supernova it presumably expands 
into a uniform ISM, ejecting roughly a Chandrasekhar mass 
Mej - I.4M0, which is characterized by the initial veloc- 
ity distribution dM^j/dv cc y^^* with k - 1. We note that, 
since SN 1006 is already in the Sedov phase, its properties are 
not sensitive to the parameter values Mgj and k of the ejecta. 
The ISM gas density po - l-4mpA^H, which is usually char- 
acterized by the hydrogen number density A^h, is an impor- 
tant parameter that strongly influences the expected SNR dy- 
namics and nonthermal emission. However, the value of A^h 
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is poorly constrained for SN 1006 from the thermal emis- 
sion observations. They yiel d a wide range of possible valu es 
from A^H = 0.05 -0.1 cm'^ towarka das & ChevaheJl"998l) to 
A^H ~ 0.3 cm-3 dPubner et a l. 2002). 

As in our early study ( Berezh ko et alJ l2002h we apply 
here nonlinear kinetic the ory of CR acceleration in SNRs 
tBerezhko et al. 1996: Berezhko & Vol k 1997) to find the op- 
timum set of physical parameters of SN 1006, mainly the 
value of A^H, which give a consistent description of the ob- 
served overall dynamics and of the nonthermal emission. As 
described in detail in these papers, that theory includes all the 
important physical factors which influence CR acceleration and 
SNR dynamics: shock modification by the CR backreaction, 
Alfven wave damping within the shock transition, selfconsis- 
tent determination of the CR spectrum, and spatial distribu- 
tion in each evolutionary phase. In addition it includes syn- 
chrotron losses of CR electrons and a determination of all non- 
thermal emission processes produced in SNRs by accelerated 
CRs (e.g. Berezhko et al. 20 02, 2003). It wa s demonstrated for 
the cases of SN 1006 (Bere zhko etalJl2002l) and Cassiopeia A 
(iBerezhko et alJ Eo03i) that the values of these key parame- 
ters (proton injection rate, electron to proton ratio and inte- 
rior magnetic field strength) which cannot be predicted theo- 
retically with the required accuracy, can be determined from a 
fit of the observed synchrotron emission data. It is important 
here that the parameter values for these SNRs, determined in 
this way, were very well confirmed by the Chandra measure- 
ments of the fine structure of the nonthermal X-rav emission 
llBerezhko et al J2003UBere7hko fc Volkl2004l). 

We take as the most reliable estimate for the distance 
d = 2.2 kpc to SN 1006 C Winkler et al, 2003) . This distance 
is roughlv 20% larger than the value c/ = 1 .8 kpc adopted in 
IBerezhko etaDll2002h . 

2.1. Ambient gas density and total explosion energy 

The SNR dynamics are illustrated in Fig. ^ In Fig. ^ we 
present the time dependence of the shock radius Rs and of the 
shock speed Vs, calculated for the optimum upstream magnetic 
field value Bo - 30 fiG (see below), and for the expected den- 
sity range of the ambient ISM. The value of the explosion en- 
ergy was taken to fit the observed size and speed at the current 
epoch f 10^ yr Since SN 1006 has already evolved into the 
Sedov phase, the explosion energy scales as Nh- The cal- 
culations correspond to £'sn/(10^' erg) = 1.9, 3.8 and 11.4, for 
A^h/(1 cm-^)= 0.05, 0.1 and 0.3, respectively. 

The scaling Esn ^ Nh is intuitively plausible, since for 
lower external density both the expansion velocity as well as 
the radius R^it) of the SNR are larger for a given explosion 
energy Esn- Asymptotically - and approximately, since there 
exists no strictly self-similar solution with particle accelera- 
tion - in the Sedov phase, we have R^it) cc (Esn/A^H)'^^?^^^- 
Therefore the ratio Esu/Nh is fixed for known t, Rsit) and 
V,(t) = dR,(t)/dt. 

The calculations shown in Fig.^ demonstrate that with the 
dependence of E^n on A^h the SNR dynamics are indeed com- 
patible with the measured values of Rs(t) and V's(0 for different 
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Fig. 1. Shock radius R^ and shock speed Vs as a function of time 
for the upstream magnetic field = 30 //G and different ISM 
number densities A^h (a), and for ISM number density A^h = 0.1 
cm"^ and different values of the upstream magnetic field ffn (b) . 
The observed size and speed of the shock llMoffettet all 19931) 
are shown as well. 



ISM densities. We note that the small differences in Ysif) for 
different A^h in Fig. la during the Sedov phase t > 400 yr are 
the result of the shock modification, which depends on A^h (see 
Fig.EJ. 

In addition, Fig.^^ demonstrates the expected weak depen- 
dence of the dynamics on the magnetic field strength Bq. For 
a higher magnetic field Bo the shock compression ratio cr is 
lower (compare Fig. |2 with Fig. lb in Berezhko et alJ ( l2002h '). 
Therefore the postshock internal pressure f 2 ~ PoK^(l - 1/c) 
is higher for lower Bo- Since in the Sedov phase the ther- 
mal energy is the dominant form of energy, the total energy is 
Esn R^Pi- Since E^^ is a constant the shock size and there- 
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Fig. 2. The overall shock compression ratio cr as a function of 
time, for three different ISM densities A^h - 0.05, 0.1, and 
0.3 cm""*, and upstream field Bo = 30 /iG. 



fore the shock speed Vs are somewhat larger for larger magnetic 
field Bo, in agreement with Fig. [2 



2.2. Amplified B-field 

In our model the interior (downstream) magnetic field is 
connected to the outer (upstream) magnetic field Bq by the 
simple relation B d = crBo, w here cr is the overall shock com- 
pression ratio (B erezhko et al. 2002) . It is also assumed that 
already the upstream field Bq is significantly amplified as a 
consequence of CR streaming, and that it therefore substan- 
tially exceeds the existing IS M value (J^ucek & BelL .2000: 
iBell & Lucekl200lUBeli2004l) . 

The shock compression ratio depends on the ambient gas 
density, as shown in Fig. |2] The behavior of the shock com- 
pression ratio cr is determined by the values of two parameters: 
the proton injection rate rj and the Alfvenic shock Mach num- 
ber Ma - Vs/ca, where Ca = Bq/ -^jAnpQ is the Alfven speed. 
The value of t he injection rate ri influenc es the critical shock 
speed V* oc T] jBerezhko & Ellisonlll999l) which separates the 
initial SNR evolutionary period, when Vs > V* and therefore 
the shock is almost unmodified, from the subsequent period 
Vs < yi of the modified shock. In our case rf - 2x 10 and 
y* «i 10^ km/s. During the period Vs < V* the shock compres- 
sion ratio goes as cr ~ \.5mI^^ ( Berezhko & Eflison 1999). 
Therefore a decrease of the ISM density A^h leads to a decrease 
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nH^^, as seen inFig.El 



of the shock compression ratio cr oc y, 

for t > 400 yr (see also Fig.^). Thus, for the case Bq = 30 fiG 
used in Fig. the downstream magnetic field strengths are 
Bd = 149, 156 and 173 fiG for A^h = 0.05, 0.1 and 0.3 cm-^ 
respectively. 



2.2.1. Spatially integrated synchrotron spectrum 

For the global determination of the effective downstream field 
Bd we compare the theoretical synchrotron spectrum, calcu- 
lated for three different values of Bq, with the observed spatially 
integrated spectrum (see Fig.|3}. 

Electrons with a power-law energy spectrum Nde) = A^e^''' 
produce a synchrotron flux Sv - Av'" with spectral index 
a = (y - l)/2 and ampHtude A oc A^B^^^^^'^ Icf. In the test- 
particle limit, y - 2 and therefore a - 0.5. Values a > 0.5, 
as observed in young SNRs, require a curved electron spec- 
trum (hardening to higher energies) as predicted by nonlin- 
ear shock acceleration models. The synchrotron emission at 
frequency v is mainly produced by electrons of energy e - 
5 V[v/(1 GHz)][(10juG)/Bd] GeV. Since the typical particle 
spectr um produced by a CR modified shock jBerezhko et alJ 
1 19961) is characterized by y > 2 at e < 1 GeV and y < 2 at 
e > 10 GeV it follows that in order to have a = 0.57, observed 
for SN 1006, numerical iteration shows that one needs efficient 
CR acceleration with a proton injection rate ?7 = 2x10 which 
leads to the required shock modifi cation, and also to a hi gh inte- 
rior magnetic field Bd> 120 /iG jBerezhko et alJ2002l) . X-ray 
synchrotron spectral measurements are required to find the op- 
timum valu e of the magnetic field strength Bd < Berezhko et alj 
I2002ll2003h : for a given fit of the synchrotron spectrum in the 
radio range the X-ray synchrotron amplitude is very sensitive 
to Bd (see Fig.Et). 

An increase of the ISM number density A^h leads to an in- 
crease of the total number of low energy CRs. Since the ra- 
dio emission is produced by low energy electrons and since for 
known SNR angular size and expansion rate the shock radius 
Rs and shock speed are proportional to distance d, we have 
Ae oc rjVsRlKgpNH ^ Tjd^Kf^pNH- Therefore the amplitude of 
the synchrotron flux A oc rjKspNud^B^J^^^^^ . Since A is fixed 
by the experiment, the electron to proton ratio changes like 
A'ep oc A^jj' for a given injection rate rj and interior magnetic 
field Bd. Therefore we have K^p = 2.6 x 10 ^ 1.1 x 10"^^ and 
4.2 X 10""* for A^H = 0.05, 0.1 and 0.3 cm""* respectively. This 
includes the renormalization factor 0.2 for the overall ion inten- 
sity. For a given A^h we have K^j, oc c/^^B^*''^''^^. Given that the 
Galactic CRs have a value K^^p ^ 10"^, the small value of Nh 
suggested by the H.E.S.S. experiment makes SN 1006 - at the 
present stage of evolution - more closely resemble an average 
source of the Galactic CRs. 

As a result the synchrotron spectrum is almost insensitive 
to A^H- The weak dependence of Sy(v) on A^h is caused by the 
different downstream magnetic field Bd - ctBq at the same Bq, 
because the compression ratio cr depends on A^h- For larger A^h 
the shock is more strongly modified. This leads to a harder CR 
spectrum at high energies e » nipC^ and therefore the spectrum 
5'v(v) at V = 10'^ - 10"^ Hz becomes harder for larger A^h, as 
can be seen in Fig.|3t. 

In addition, an increase of the magnetic fie ld Bh leads to a 
decrease of the electron energy e\ oc r^B^^ jBerezhko et alJ 
2002), where the electron spectrum undergoes a brake: for 
e > e\ synchrotron losses in the downstream region are sig- 
nificant and lead to a steep overall (spatially integrated) elec- 
tron spectrum A^e e"''"\ where j is the power law index 
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log V, Hz 

Fig. 3. Synchrotron emission flux as a function of frequency for 
the same cases as in Fig.^ The observed X-ray ("Hamihon et alj 
[l986; Allen et al. 1999.) and radio emission (■Revnolds..l996j) 
fluxes are shown. 

of the CR spectrum at the shock front. Therefore the syn- 
chrotron spectrum has a break at the corresponding frequency 
v\ oc ej^Z?d oc B^^ , which decreases with increasing Bd, as one 
can see in Fig.|3] 

An increase of the magnetic field requires a reduction of 
the amplitude of the electron momentum distribution since the 
radio electrons are not subject to synchrotron cooling at the 
present epoch. This leads to a decreased electron/proton ra- 
tio /Tep oc Bq" ' in the energy range where radiative cool- 
ing is unimportant. Most importantly however, the cooling re- 
gion of the synchrotron emission spectrum v > v\ will be 
lowered in amplitude by this decreasing electron/proton ratio. 
Comparison with the experimental X-ray data shows that the 
optimum magnetic field value is about Bq - 30 [iQ, and cor- 
responds to a downstream field k. 150 juG. This is some- 
what larger than the B^ ^ 120 AfG quite conservatively esti- 
mated earlier (Berezhko et alJ2002ll2003h . and is in agreement 



with the field amplification that is implied by th e filamentary 
structures in hard X-ray s, observed with Chandra jBamba et alJ 
l200:4IVolk etal .'2005') (see also below). 

As it was already noted (Berezhko et al. 2002), it is only 
the highest energy part of the X-ray spectrum corresponding to 

V > 10'^ Hz that is of nonthermal origin. At lower frequencies 

V <K 10'^ Hz the spectrum is dominated by thermal emission. 
Therefore the theoretical spectrum for v <K 10'^ Hz should be 
substantially below the experimental values, as it is in Fig. 3. 

2.2.2. Overall synchrotron morphology 

We now compare this globally determined magnetic field value 
B() - 30 juG with the local Chandra results by plotting the full 
numerical solutions for the overall SNR morphology, as it re- 
sults from selfconsistent time-dependent kinetic theory, jointly 
with the sharpest of the experime ntal brightness p rofiles ob- 
tained by the Chandra observers ( Bamba et al. 2003) in Fig.|^, 
where the values of I do f characterizing the quality of fit are 
also shown. Since this profile obviously lies on a high superim- 
posed background, we have subtracted from the data the con- 
stant brightness value which corresponds to the far upstream 
region. Note that the initial Chandra data were presented as a 
function of the angular radial distance ( Bamb a et alJl2003l) . 
The transformation to the linear projection distance p - i/zd 
has been achieved by assuming a distance of li = 2.2 kpc to 
SN 1006 (see above). We use in Fig.l^lSlAi// instead of ip, tak- 
ing the initial point Ai/r = at the theoretically predicted shock 
position. The agreement of global morphology and local pro- 
files is reasonable for all considered ISM densities. 

Since the absolute values of the experimental brightness 
profiles and the actual shock position are not known their val- 
ues were used as adjusting parameters in the fitting procedure. 

Note that the radial profiles shown in Fig. |^ depend on 
A^H as a result of the dependence of the downstream magnetic 
field Bd - ctBq on A^h: for higher A^h the compression ratio 
cr is larger (see Fig. 2) and leads to a higher Bd and therefore 
to a thinner profile. In addition, the radial profile depends on 
the downstream plasma speed U2 - Vsjcr (Berezhko & Volk 
2004). The lower value of ui, that is associated with a larger 
A^H, makes the profile thinner 

As can be seen from Fig.|^, the magnetic field value Bq - 
20 slightly overestimates the profile width. Together with 
the spectral fit (Fig. 13)) we conclude that Bq - 30 //G (Bd ~ 
150 yuG) is a realistic estimate. 

2.2.3. Local fit of the synchrotron profile 

X-ray brightness profiles can also be used for an independent 
determination of the interior magnetic field value that does not 
require the nonlinear theory (Fig.|5}. 

It has been shown jBerezhko et alJ2 003': Ber ezhko & Volkl 

I2OO4I) that in the case of strong synchrotron losses, which 
in fact occur for the electrons responsible for the X-ray syn- 
chrotron radiation, the downstream diffusion length over a syn- 
chrotron loss time mainly determines the three-dimensional 
thickness of the X-ray emission profile. The effective magnetic 
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Fig. 4. Projected radial dependence of the non-thermal X-ray 
brightness in the energy range 2 to 10 keV, calculated for the 
same cases as in Fig. ^ together with the Chandra da ta corre- 
sponding to the sharpest profile in jBamba et alJ2003l) . 



field can then be expressed in terms of the width li of the radial 
synchrotron emissivity profile qv{r) according to the simple re- 
lation 
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where m^, e and ro denote the electron mass, charge, and clas- 
sical radius, respectively. The actual value of I2 is determined 
from the observed width of the projected radial brightness pro- 
file 
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where p is the projected distance from the remnant center and 
the integration is performed along the line of sight. 




p/Rs 

Fig. 5. The same Chandra 2-10 keV radial continuum bright- 
ness profile as in Fig.|3] fitted to the projection of the exponen- 
tial emissivity profile described by the model of Eq. (|2jl. The fit 
has ax^/dof = 1.2. The dashed lines indicate the Icr deviation. 



Since the synchrotron emis sivity can be repre sented in the 
form qvir) — q2 exp[(r-7?s)//2] ( iBerezhko et al.l2002) . the pro- 
jected profile is expressed in the analytical form 



Jy 



qi exp(-f ) erfi (f), 



(2) 



where t - -\/R^-fP- /(p ^|2RJ2), and where erfi(x) is an 
imaginary error function. Note that the profile Jvip) can be 
expre ssed in terms of elementary functions (Berezhko & Volk 
I2OO4I) . if we use the expansion of erfi(f) near the point f = 0. 

As shown in Fig. |5] the experimental Chandra values are 
quite well fitted by the analytical profile, described by Eq. 
with an emissivity width h - 13.7(-i-8.2 - 6.2) x 10"^ pc, that 
results in the effective magnetic field - 191(+82 - 62) juG. 
This range of downstream B-field strengths includes the values 
following for the three densities considered in Fig. on the 
basis of an upstream field of Bo = 30 yuG which we shall use 
for the following. 

Taking into account these results together with the consid- 
eration presented in Fig. 3 we conclude that Bq - 30 //G and 
Bd ~ 150 fiG are the most appropriate values consistent with 
the existing data. 

3. Nonthermal energy density 

As we have seen, fisn °^ Nh in the Sedov phase, for given SNR 
size and expansion rate. Because these parameters are indepen- 
dent of time, it is necessary to choose this proportionality gen- 
erally. Fig.^demonstrates that this approach is roughly consis- 
tent for all phases. The resulting values for in Fig.|6lare also 
calculated on this assumption. 

In Fig.|5]we present the relative CR energy content E^/Esn 
as a function of time, calculated for So = 30 fiG and three dif- 
ferent ISM number densities A^h- It increases with time, reaches 
a peak value at the epoch f ^ 3000 yr and then slowly decreases 
due to an increasing dominance of adiabatic cooling. In the 
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Fig. 6. Total energy EdE^^ of accelerated particles in the SNR, 
normalized to the total hydrodynamic explosion energy fisn, as 
a function of time for the different ISM densities A^h shown in 
Figs. 1-5, assuming spherical symme try. E^ is renormalized by 
a factor of 0.2 (cf. IVoUc et aTll2003h . The upstream magnetic 
field strength is Bq = 30 ju G. 



early free expansion phase f ^ 100 yr, when nonlinear shock 
modification is not essential (see Fig.|2ji, the shock produces a 
power law CR spectrum / oc MnjPinjP where A^i„j = rjNu is 
the number of gas particles injected into the acceleration from 
each unit of volume intersecting the shock front and pi„j oc 
is their momentum. Since the strong unmodified shock always 
produces a CR spectrum of the same shape, the CR energy con- 
tent Ec oc Nc at this stage grows proportionally to the total 
number A^c of accelerated CRs if we neglect weakly time de- 
pendent factor ln(/?tnax/'M/>c), where /^max is maximum CR mo- 
mentum. Assuming spherical symmetry for the moment, and 
starting from the obvious relation 

dNc oc NuRlvldt 



and ta king into account the shock expansion law (eg. lChevaheJ 
Il982h 



(3) 



R, oc {EljNn)'"t'l' 
we have 



E^ocE^Nll'f^^ 



Then we have a residual density dependence at this early phase 
Ec/Esn oc N^'' if we use £511 Nh- 

In the transition to the Sedov phase (which happens earlier 
for larger A^h), Ec/Esd slowly reaches a maximum as a function 
of time which is roughly independent of A^h- This is the result 
of the nonlinear limitation by the total amount of energy E^n- 
During the Sedov phase the shock speed ¥^(1) depends only on 
E^n/Nii, as we have seen before (see also Fig.[^), whereas the 
Alfvenic shock Mach number Ma - Vs/ca ^ ^/Nh increases 
with A'h, because the Alfven speed oc 1/ -\/Nii- This leads 




Fig. 7. Integral jr^-decay (thick lines) and IC (thin lines) y-ray 
fluxes from the NE half of the remnant as a function of y-ray 
energy for Bo = 30/iG, and for the same different densities as 
in Fig. n The CANGAROO-I (open circles with error bars) 
and CANGAROO-II (solid circle with error bars) flux from 
the NE rim (Hara 2002) and the cor respondingJI.E.S.S. up per 
hmit for the CANGAROO position jAharonian et alJ2005l) are 
shown as well. 



to a stronger shock modification (see Fig.|5J, that in turn pro- 
vides a more efficient CR production Ec/Esn for larger A^h- In 
spherical symmetry E^ can reach more than 50 % of E^n in this 
way. In reality, however, this spherical symmetry does not ex- 
ist, because for a SN explosion in a uniform medium, as for 
SN 1006, the external magnetic field will also be uniform. This 
leads to a severe reduction of the injection rate around those 
parts of the shock surface that are roughly parallel to the exter- 
nal field and requires a renormalization of th e overall accelera - 
tion efficiency by a factor of the order of 0.2 ( Volk et al. 20031). 
This dipolar picture of the synchrotron emission has been con- 
firr ned in an elegant analysi s of recent XMM data for SN 1006 
by 'Rothenflug et alJ ( l200 4l). Altogether this leads to a rough 
proportionality E^ oc Nh until particles start to leave the rem- 
nant. As one can see from Fig. |6j the normalization of Ec by 
Esn removes the main density dependence: the relative differ- 
ence £'c(A^H = 0.3)/£'c(A^H = 0.05) - 1 of the values equals 
about 1 1 , whereas the relative difference of the normalized val- 
ues Ec/Esn is about unity. 



4. Gamma-ray fluxes 

The previous considerations determine the values of the proton 
injection rate, the electron/proton ratio and the magnetic field 
strength Bq. The only important parameter which cannot be de- 
termined from the analysis of the synchrotron emission data is 
the external density A^h- Therefore we now compare calcula- 
tions for different A^h with the existing y-ray measurements. 
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Figs.0and|S]show the calculated integral y-ray fluxes from 
the North-Eastern (NE) half of the remnant, in order to compare 
with the reported CANGAROO emission from the NE rim and 
the corresponding upper limit from H.E.S.S. These calculated 
fluxes correspond to 50 percent of the entire theoretical y-ray 
flux from SN 1006 under our assumption of dipolar symmetry 
(Volketal. 2003). 

The key point is that a decrease of the ISM density leads to 
a considerable reduction of the expected /r^-decay y-rays (see 
Fig.0: in the Sedov phase the ;7r*^-decay y-ray flux F'^ oc A^h^^c 
is proportional to the gas density and to the CR energy content 
Ec- This roughly implies oc A^^, as shown in the previous 
section and in agreement with the numerical results. In addi- 
tion, an increase of the gas density leads to an increase of the 
Alfvenic shock Mach number and of the shock modification 
due to CR backreaction. This implies in turn a harder CR spec- 
trum and correspondingly an increase of the y-ray flux at the 
highest energies. Only this last factor influences the IC y-ray 
flux. Therefore it is only weakly sensitive to the gas density as 
Fig.EJalso indicates. 

To allow for a simple reading of the figures, we present the 
same data in two separate figures, Fig.0and Fig.|8t, showing 
in Fig.0the theoretical 7r"-decay F" and IC F^ integral fluxes, 
whereas in Fig.|8t we show total fluxes Fy - Fy + F^^ and Fjp. 

In Fig. 1^1 (see online version of this paper) we show the 
total (TT^-decay + IC) and the IC integral y-ray fluxes, again 
from the NE half of the remnant as a function of y-ray energy 
for the same cases as in Fig.^ Also the total y-ray flux varies 
strongly with the ISM density (Fig0. 

A modification of the magnetic field strength would not 
change much the 7r°-decay y-ray flux, and thus the total y-ray 
flux, whereas the IC flux decreases considerably with increas- 
ing magnetic field strength. Fig. [SJi illustrates this fact, even 
though we consider the value Bo ~ 30//G as experimentally 
quite well determined. 

Note that an Inverse Compton y-ray emission scenario in a 
low magnetic field of order B^/ = 4± 1 /zG dTanimori et alJ20oTh 
is clearly excluded by the H.E.S.S. upper Hmits, and this is fully 
consistent with the measured and expected field amplification. 

T he non-detection of SN 1006 by the H.E.S.S. experi- 
ment jAharonian et alJl2005h implies that the hydrogen den- 
sity A^H is lower than 0.1 cm~^, and excludes the high value 
A^H = 0.3 cm"^ which happened to fit the CANGAROO data. 
As clearly shown in Fig.|8t, at 1 TeV and for Bq - 30yL(G the 
H.E.S.S. upper limit is still about three times larger than the 
total y-ray flux from SN 1006, expected for A^h = 0.05 cm"^, 
which is at the lower end of the range 0.05 < A^h ^ 0.3 cm"^ 
of plausible ambient densities. 

Since the IC y-ray emission is quite insensitive to the ISM 
density, the lower limit for the expected y-ray flux at TeV en- 
ergies - as derived from the integrated synchrotron flux and the 
field amplification alone - is about a factor of five lower than 
the H.E.S.S. upper limit. Such a low y-ray flux from SN 1006 
would be expected if the ambient ISM number density A^h was 
even considerably lower than 0.05 cm"^. 



5. Summary 

Our reexamination of the most relevant set of physical and as- 
tronomical parameters of SN 1006 led us to chose the some- 
what greater source distance of 2.2 kpc, determined by more 
recent measurements, and to establish the interior magnetic 
field in SN 1006 more precisely at » 150 juG. Most im- 
portant, however, is another fact: the lack of a TeV signal from 
SN 1006 that follows from the non-detection by the H.E.S.S. 
instrument does not invalidate the theoretical picture which 
gives a consistent description of the nonthermal emission char- 
acteristics; it rather implies a constraint on the ambient gas 
density A^e < 0.1 cm"^. This point was discussed in detail. 
The renormalization of the y-ray flux on which this result is 
partly based has been recently confirmed experimentally. The 
well-defined morphology makes this remnant the simplest and 
best understood of its kind. In conjunction with the amplifi- 
cation of the magnetic field by the accelerated cosmic rays, 
the H.E.S.S. non-detection completely rules out the low-field 
Inverse Compton scenarios that were popular until a few years 
ago. 

In the case of SN 1006 and in other similar cases y-ray 
measurements will at the same time give a reliable estimate of 
the ambient ISM density. 

Quantitatively (see Fig.0and Fig.|8^), the present H.E.S.S. 
upper limit is only a factor of about three above the total TeV 
flux expected for A^h = 0.05 cm"^, which is at the lower end 
of the range 0.05 < A^h ^ 0.3 cm"^ discussed in the literature. 
If the actual ISM is so diluted that A^h < 0.05 cm"^ then the 
lowest expected TeV y-ray flux is by a factor of five below the 
present H.E.S.S. upper limit. 

Given the arguments above, it is certainly worth attempting 
the detection of this important source in a deep observation of 
about 200h with the H.E.S.S. experiment. 
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, TeV 

Fig. 8. (On-line-only) Total (7r°-decay + IC) (thick lines) and 
IC (thin lines) integral y-ray fluxes from the NE half of the 
remnant as a function of y-ray energy for the same cases as 
in Fig. n The CANGAROO-I [open circles with error bars) 
and CANG AROO -II (solid circle with error bars) flux from 
the NE rim llHarai2 002) and the cor responding H.E.S.S. u pper 
limit for the CANGAROO position dAharonian et alJ2005l) are 
shown as well. 



